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Abstract 


In this paper the line-spring model developed by Rice and Levy is 
used to obtain an approximate solution for a cylindrical shell contain- 
ing a part-through surface crack. It is assumed that the shell con- 
tains a circumferential or axial semi-elliptic internal or external 
surface crack and is subjected to a uniform membrane loading or a uni- 
form bending moment away from the crack region. To formulate the shell 
problem, a Reissner type theory is used in order to account for the 
effects of the transverse shear deformations. The stress intensity 
factor at the deepest penetration point of the crack is tabulated for 
bending and membrane loading by varying three dimensionless length para- 
meters of the problem formed from the shell radius, the shell thickness, 
the crack length, and the crack depth. The upper bounds of the stress 
intensity factors are provided by the results of the elasticity solution 
obtained from the axisymmetric crack problem for the circumferential 
crack, and that found from the plane strain problem for a circular 
ring having a radial crack for the axial crack. Qualitatively the 
line-spring model gives the expected results in comparison with the 
elasticity solutions. The results also compare well with the existing 
finite element solution of the pressurized cylinder containing an inter- 
nal semi-elliptic surface crack. 


1. Introduction 

In recent years there has been some renewed interest in the line- 
spring model which was developed in [1] for obtaining an approximate 
solution of a plate containing a part-through surface crack. There are 
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a number of reasons for this. First, the accuracy of the results 
obtained from the model turned out to be better than that shown by the 
early comparisons with the solutions found from the finite element and 
the alternating methods [2-6] (see, for example, [7]). Secondly, the 
technique appears to have the potential for important applications to 
a great variety of shell structures of rather complex geometries with 
a relatively small computational effort. Finally, it can be quite 
useful to study certain aspects of the part- through crack problem in 
the presence of large scale plastic deformations (see, for example, 
the interesting recent work by Parks [7,9], and [8] and [10]). 

In this paper the elastic problem for a relatively thin-walled 
cylinder containing a semi-elliptic part-through crack is considered. 

It is assumed that the crack lies in a plane perpendicular to or con- 
taining the axis of the cylinder and may be an external or an internal 
surface crack. In formulating the problem, the cylinder is approxima- 
ted by a shallow shell and the effect of transverse shear deformations 
are taken into account [11,12]. The edge-cracked strip results used 
in the line- spring model are obtained from an integral equation solu- 
tion given in [13]. 

The stress intensity factor for a part-through axial crack located 
inside the cylinder is given in [14-16] where in [14] and [15] the 
finite element and in [16] the boundary integral equation method is 
used to solve the problem. The results found in this paper are com- 
pared with those given in [14] as well as the related plane strain and 
axisymmetric elasticity solutions. The stress intensity factors 
obtained from the elasticity solutions for a ring with a radial crack 
under plane strain conditions and for a cylinder containing an axi- 
syiranetric circumferential surface crack provide upper bounds for the 
results corresponding to an axial and a circumferential surface crack 
of finite length. 
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2. General Formulation 






The part-through crack geometry for the cylindrical shell under 
consideration is shown in Figure 1. It is assumed that the external loads 
are symmetric with respect to the plane of the crack. Thus,, the only 
nonzero net ligament stress and moment resultants which have a con- 
straining effect on the crack surface displacements would be the membrane 
resultant N-j^ and the moment resultant M-j-j. The basic idea underlying 
the line-spring model consists of (a) representing the net ligament 
stresses by a membrane load N and a bending moment M, and the crack sur- 
face displacements by a crack opening 6 and a relative rotation e, all 
referred to the midplane of the shell and continuously distributed along 
the length of the crack, (b) approximating the relationship between 
(N,M) and ( 5 , 0 ) by the corresponding plane strain results obtained from 
the solution of an edge-cracked strip or a ring, and (c) reducing the 
problem to a pair of integral equations for the unknown functions N 
and M or 6 and 0 by using the boundary and the continuity conditions 
for the shell in the plane of the crack. 

In the formulation of the crack problem for the shell, the deriva- 
tives of the crack surface displacement and rotation are used as the 
unknown functions which are defined by 


^ u(+0,y) = G^Cy) , ■- = G 2 (y) 


(la,b) 


The notation and the dimensionless quantities used in the formulation 
are given in Figure 1 and in Appendix A. It is shown in [17] that the 
general problem for a symmetrically loaded shell containing a through 
crack may be reduced to the following system of integral equations: 


1 

-1 


G-] (t) 


dt + 


1 

» 

-1 


2 k^ .(y,t)G.(t)dt = 2TTF^(y), -1 < y < 1, 


(2a) 


l-v2 


-1 


G2(t) 


dt + 


-1 


1 a -l<y<l 


(2b) 
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CO 


k-|i(yjt) - 


[2 Z ct2Q.(a) - 1] sina(t-y)da , 
1 




2a^ z N -(a)sina(t-y)da , 
1 


i<2i(y>t) = 


^22 


2 f 

TT 


4 p5(.m^-va2)Q.(oi) 






4 p5(mi-\'a^)N;:(a) 


- ic(l-v)2ar-j + (l-v2)/2]sincx(t-y)da 


(3a-d) 


subject to 
1 


G-|(y)dy = 0 


G2(y)dy = 0 


-1 


-1 


(4a, b) 


The problem is formulated as a stress disturbance problem in which a 
homogeneous stress solution for the uncracked shell is separated through 
a superposition and it is assumed that the stress and moment resultants 
applied to the crack surfaces are the only external loads. Thus, F-j 
and p 2 appearing in (2) are 


Fi(y) = -l<y<l . (5a, b) 


The parameters r-,, m., and p., (j=l,...,4) are functions of the trans- 

I o 

form variable a and are given by 



2 

k ( I - v )'' 


( 6 ) 


= Pj + a2 , Re(mj) < 0, (j=l 


,4) 


(7) 
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and pp...,p^ are the roots of 


- icA^P^ + {ZKX^X^a^ - 2^^201^ + + (2KX^X|a^ - KX'^a^ 

- KX!|a2 + 2X^ - 2x2x|)a2p + (x|-X^) o'* = 0 (8) 

From (6-8) it may be shown that for large values of |a| we have 

r,(») = -|c.| [1 + + 0(^)] , (9) 

2 

nij(a) = -|al [1 + 2^ - 8^ + 

where the roots p. of the characteristic equation (8) are bounded for 

V 

all values of a. 

The functions Q- and N-, (j=l,...,4) which appear in the kernels 

J u 

(3) are found from 

Rj(a) = i[Qj(a)f-|(a) + Nj(a)f2(a)] , (j = 1,...,4) , (11) 


where 


^k<“> = 


1 

* 

-1 


G,^(t)e^'“^ dt, (k=l,2) 


and Ri,...,R^ are obtained from 

4 

E m-R-(a) = 0 j 
1 J J 


X^mj-X^a^ 


Rj(a)pJn,i 

_ i(l-v)K 

(KPj-l)(X^m4-X^a^) 

2aX^ 

R.(»)pjn,j 


((cPj-l)(x|m^-x^a'^) 

- a 

X^ “ 


( 12 ) 


(13a-d) 
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The formulation given above refer to a shallow shell containing 
a crack along the principal plane of curvature coinciding with X2X2 
plane (Figure 1). The principal radii of curvature and R2 are 
defined by 


1 32Z 1 _ 3^2 

R^ ~ " ’ R2 ■ " IX^ ’ 


(14a, b) 


where Z(XpX2) is the distance of the point on the middle surface to 
the tangent plane X-jX2. Thus, for the circumferential crack shown in 
Figure la, R2=R and R-|=o> (giving Xi=0), and for the axial crack, shown 
in Figure lb R.j=R and R2=” (giving X2=0). 

Let now 


Nil = . Mil = (15a,b) 

be the uniform membrane load and the bending moment applied to the 
shell away from the crack region and N(X2) and M(X2) the stress and 
moment resultants which are equivalent to the net ligament stresses in 
-a<X2<a or -l<y<l. The "input" functions which appear in the integral 
equations (2) may then be expressed as 

Fl(y) = (“CTco'*''^) > ^ (-nico+m) , (16a, b) 


for the crack located on the outside and 


Fi(y) ^ (“CT(„+cj) , F2(y) g£ (m„,“m) , 
for the crack located inside the cylinder^ ' where 


_ 6M00 


CToo = > nioo = -j^ 


(17a, b) 


(18a, b) 


TIT] 

' ' In both cases the applied moment M is such that the crack tends 
to open under bending of the shell, and the net ligament moment 
is assumed to constrain the crack surface rotation, hence the 
change in sign of F2 in (16) and (17). 
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and 





N(av) 
h ~ h 


eMCXp) 
, m(y) = ^ 


_ 6M(ay) 
h2 * 


09a, b) 


The stresses a and m are linearly related to the crack surface dis- 
placement a u(+0,y) = 5/2 and rotation 3x(+0>y) = 0/2. This relation- 
ship may be obtained from the related plane strain problem by expressing 
the rate of change of the potential energy in terms of the crack closure 
energy and the change in gross compliance as follows: 


l-v2 

E 


iv2 _ 1 r t 3 5 ■ 
K -^[ah^ + 


mh^ 39i 

6 at-* 


C20) 


where K is the total mode I stress intensity factor at the crack tip 
and L is the length of the edge crack. If we now let 

K = /h (a g^ + m gj^) , (21) 


from (20) we obtain 

a(y) = E[Y^^(y)u(+0,y) + ’ 

m(y) = 6E[Y(j^(y)u(+0,y) + ’ (22a, b) 

where + and - signs are to be used for the outer and the inner cracks, 
respectively and 

a “bb _ 1 “tt 

^tt " h(l-v^) ~ ’ ^bb " 36(1 -v2) a ’ 

1 “tb _ a “bt 

^tb ■ ■ 6(1 -v'2) “a~ ’ ^bt ■ ■ 6h(l-v^) A ’ 

^ " "‘tt “bb " “tb ’ • (23a-e) 

I L 

“ij "" F-'o ^i^j 
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The crack depth L is assumed to be a known function of y (Figure 1). 
Referring to the definitions (1), u and 3^^ may be expressed as 


r 

u(+ 0 jy) = J 6 -|(t)dt , 3 j^(+ 0 >y) = G2(t)dt 
-1 -1 


(25a, b) 


Substituting from (25) and (22) into (2) the final form of the integral 
equations is found to be 

/ / T r' G,(t) 

J 


-1 


J_ 

2ir 


-1 

.y 


[k-ji(y,t)G^(t) + ^ ~ ~T ’ » 

^ G,(t) 


;Ybt(y)f Gi(t)dt - Ybf^(y)| G2(t)dt + I -|iy-clt 

-1 -1 -1 


2irh 


m„ 


Cl<2i(y»'t)Gi(t) + k22(y,t)G2(t)]dt = + ^ > -l<y<l 


-1 


(26a, b) 


where the upper (i.e. , -) and the lower (i.e., +) signs are to be used 
for the outer and the inner crack, respectively. 


3. Compliance Coefficients 

The functions g^ and g^^ which appear in (21) and which give the 
membrane and bending components of the stress intensity factor are 
obtained from the corresponding plane strain crack geometry. 

For the circumferential crack, the appropriate geometry is that 
of an infinite strip with an edge crack. On the other hand, for the 
axial crack the proper plane strain problem would be that of a ring 
having a radial edge crack. In a recent study the ring problem was 
formulated in terms of a singular integral equation [18]. The results 


I* 
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given in [18] show that for cylinders with values of h/R which may be 
considered a "shallow shell", the ring results are reasonably close to 
the strip results. Also for small values of h/R the convergence of 
the numerical solution of the ring problem is not very good. Hence, 
the complete parametrization of the problem for the purpose of obtain- 
ing g^ and g^^ (which would be functions of h/R as well as L/h) becomes 
rather complicated. In this paper, therefore, the edge-cracked strip 
results will be used for both the axial and the circumferential crack 
problem. 

For the strip the functions g^ and gj^ are obtained from the results 
given in [13] which are valid for 0<L/h<0.8 and may be expressed as 

g^(c) = (1,1216 + 6.5200^2 - 12.3877?'^ + 89.0554?6 

- 188.6080^9 + 207.38705^9 - 32.0524512] , 


gb(?) = (1.1202 - 1.88725 + 18.014352 - 87.38515^ 

+ 241.91245'^ - 319.94025^ + 168.010559) (27a, b) 


where 5 = l(X 2 )/h = L(ay)/h. From (27) and (24) the functions a^-j, 
(i,j=t,b) may be determined as follows: 


‘tt 


5^ S C'I> . 
n=0 


‘‘bb 


2 5" 


“tb ■ “bt 


= 


19 (n) n 

n=0 


(28a-c) 


The coefficients 


c!'!) 

IJ 


are given in Table 1. 
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4. Solution for the Cylindrical Shell, 

The solution of the problem is obtained for a uniform membrane 
loading N„ and for a bending moment Moo applied to the shell away from 
the crack region and for the Poisson's ratio v = 0.3. Even though 
L(X2) = L(ay) describing the crack shape can be any single-valued func- 
tion, the problem is solved only for a semi-elliptic surface crack 
given by 

L = Ljj/RyiTF = . (29) 

The solution of the integral equations (26) is of the form 

(t) 

G,-(t) = -- - A , , (i=l,2) C-l<t<l) (30) 

^ (l-t^)^ 

where (j>^ and <^2 bounded functions. The functions may be deter- 
mined from (26) to any desired degree of accuracy by using the Gauss- 
Chebyshev integration procedure [19]. After obtaining ij)-j and 412 the 
unknowns a and m representing the net ligament stresses may be deter- 
mined from (22) by using (23-25), (27) and (28). The stress intensity 
factor K(y) may then be obtained from (21) and (27). For a Poisson's 
ratio V = 0.3 and for various crack geometries and loading conditions 
the calculated results are shown in Figures 2-7 and Tables 2-11. 

Tables 2-9 give the normalized stress intensity factor at the deepest 
penetration point y=0, 1=1^ of a semi-elliptic surface crack in a 
cylindrical shell under uniform membrane loading or bending. The 
normalizing stress intensity factor k is the corresponding value 
for the plane strain problem under tension or bending and is given 
by 

Kn K Ln 

= . ( 31 ) 

VTT vir h 

for membrane loading, and 
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(.32) 


k 


0 



6M^ 


^0 



for bending. 

Figures 2 and 3 show the comparison of the shell results with the 
stress intensity factors obtained from the corresponding axisymmetric 
and plane strain problems. As CR-j/Rq) ^ 1 the shell results approach 
the flat plate solution kp [21] having a part-through semi -elliptic 
crack of the same geometry and relative dimensions. It may be noted 
that, as expected, the shell stress intensity factors are generally 
smaller than the corresponding two-dimensional values. Even though 
the shell results are given for 0.74<(R-{/Ro)<l> because of the nature 
of the theory used in the shell analysis, namely the shallow shell 
theory, for (R^/Rg) < 0.9 the results may not be very accurate. From 
Figure 3 one may also observe that the difference between shell and 
the plane strain results decrease with decreasing crack depth (Lg/h)*. 

Some sample results for the distribution of the stress intensity 
factor along the crack front are given in Tables 10 and 11 and in 
Figures 4-6. The normalization factors kg used in these tables and 
figures are also those given by (31) and (32). The variable 4 > used 
in the presentation of these results is the usual parametric angle of 
the ellipse shown in the insert of Figure 4. For small values of ({> 
the stress intensity factors given by the line-spring model are neither 
reliable nor meaningful and therefore are not presented. 

The only shell results which exist in literature and which are 
obtained by using a method other than that of the line-spring are the 
stress intensity factors in a pressurized cylinder containing an 
internal semi-elliptic axial crack [14-16]. Figure 7 shows the com- 
parison of the stress intensity factors obtained from the line-spring 
model and those given in [14] which are found by using the finite 
element method. The parameter (j> is again defined by the insert in 

' The plane strain cylinder results given in Figure 3 are obtained 
from [18] and the axisymmetric crack results shown in Figure 2 
are from [20]. 
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Figure 4. The stress intensity factor ratio F shown in Figure 7 is 
def i ned by 

■ pRj ’ 

where K = k/n is the stress intensity factor along the crack front, p 
is the internal pressure and Q = [E(k)] , E being the complete elliptic 
integral of the second kind. The results given in Figure 7 include 
the effect of the pressure p acting on the crack surface. Considering 
the gross approximations involved in the formulation of the problem by 
using the line-spring model, and the fact that the finite element 
results themselves may contain a few percent error, the agreement between 
the two results seems to be quite good. The plane strain results 
given in Figure 3 suggest that the accuracy of the results given by 
the line-spring model could perhaps be improved further if the ring 
rather than the flat plate solution is used to derive the functions 
9t ^b express the stress intensity factor (see equations (21) 
and (27)). 
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APPENDIX A 


The notation and dimensionless quantities (Fig. 1) 


X = 


T ‘ y 


z = 


u = 


_ 


= ^ 


W 

, w = 


“ 3-| » 3y ~ 3 


2 ’ 


N-ji Npo N-|2 

N = N = -- - N = -■ 

*\x hE ’ ‘Vy hE ’ ‘V hE 

M-ji M22 M-j 2 

^xx = FE ’ " Fe ’ % = FT 

V] V 2 

v = — V = — 

\ hB ’V hB ’ 


a; = i2(W) = i2o-v2) ^ 


® ° Ijrfcr > = i2(i-v2) ^ ^ 


Up U 2 , W: components of the displacement vector, 

3p 32 : rotations of the normal, 

N--, (ijj=l,2): Membrane stress resultants 

* J 

M--, (i,j=l,2): Moment resultants 

* vl 

Vp (i=l,2): Transverse shear resultants 

Rp R 2 : Principal radii of curvature 


(A.l) 

(A.2) 

(A. 3) 
(A.4) 

(A. 5) 

(A.6) 

(A.7) 
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(n) 

Table 1. Coefficients C. . which appear in eqs. (28) 

0 


n 

r(n) 

^tt 

'^tb 

c(n) 

^bb 

0 

1.9761 

1.9735 

1.9710 

1. 

11.4870 

-2.2166 

-4.4277 

2 

7.7086 

21.6051 

34.4952 

3 

15.0143 

-69.3133 

-165.7321 

4 

280.1207 

196.3000 

626.3926 

5 

-1099.7200 

-406.2608 

-2144.4651 

6 

3418.9795 

644.9350 

7043.4169 

7 

-7686.9237 

-408.9569 

-19003.2199 

8 

12794.1279 

-159.6927 

37853.3028 

9 

-13185.0403 

-988.9879 

-52595.4681 

10 

7868.2682 

4266.5487 

48079.2948 

11 

-1740.2463 

-2997.1408 

-25980.1559 

12 • 

124.1360 

-6050.7849' 

6334.2425 

13 


8855.3615 


14 


3515.4345 


15 


-11744.1116 


16 


4727.9784 


17 


1695.6087 


18 


-845.8958 
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Table 2 . Normalized stress intensity factor k/ko at the deepest 

penetration point L=Lq, y=0 of an outer semi elliptic cir- 
cumferential crack in a cylinder under uniform membrane 
loading Noo; = [12(l-v2)]^'+ a//^, v=0.3. 




^0 - 

u./in 



■-o - 

u.^tn 


^2 


a=2h 

a=4h 

a=8h 

■ilBI 

a=2h 

a=4h 

a=8h 

0 

0.817 

m 

0.930 

0.961 


Bi 

0.741 

0.837 

0.5 

0.817 

IH 

0.930 

0.961 


tjm 

0.742 

0.837 

0.75 

0.816 

0.882 

0.930 

0.961 



0.742 

0.837 

1.0 


0.880 

0.929 

0.960 


0.626 

0.741 

0.836 

1.5 


0.876 

0.926 

0.959 


0.620 

0.736 

0.833 

2.0 



0.922 

0.956 



0.727 

0.827 

4.0 



0.893 

0.939 



0.670 

0.784 

6.0 




0.916 




0.728 

8.0 




0.893 




0.676 



Lo = 

0.6h 



1-0 = 

0.8h 


0 

0.245 


0.451 

0.582 

0.073 

■mi 


0.216 

0.5 

0.248 

IM 

0.454 

0.583 

0.074 


Mg 

0.219 

0.75 

0.250 

0.341 

0.455 

0.585 

0.076 


|M 

0.220 

1.0 


0.341 

0.455 

0.585 


0.109 

0.154 

0.221 

1.5 


0.341 

0.453 

0.583 


0.112 

0.157 

0.223 

2.0 



0.448 

0.577 



0.158 

0.224 

4.0 



0.408 

0.532 



0.158 

0.214 

6.0 




0.476 




0.197 

8.0 




0.428 




0.182 
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Table 3. Normalized stress intensity factor k/ko at the deepest 
penetration point L=Lq, y=0 of an outer semi -elliptic 
circumferential crack in a cylindrical shell under uniform 
bending moment Mo,. 




"o - 

u.zn 



■-0 = 

u.^fn 


h 

a=h 

a=2h 

a=4h 

a=8h 

a=h 

a=2h 

a=4h 

a=8h 

0 

0.804 

0.875 

0.926 

0.959 

0.441 

0.579 

0.710 

0.819 

0.5 

0.804 

0.875 

0.926 

0.959 

0.443 

0.581 

0.712 

0.819 

0.75 

0.803 

0.874 

0.925 

0.958 

0.443 

0.580 

0.711 

0.819 

1.0 


0.872 

0.924 

0.958 


0.578 

0.709 

0.818 

1.5 


0.867 

0.921 

0.956 


0.570 

0.703 

0.814 

2.0 



0.916 

0.953 



0.692 

0.806 

4.0 



0.884 

0.934 



0.621 

0.753 

6.0 




0.909 




0.686 

8.0 




0.883 




0.624 



'-0 = 

0.6h 



'-0 = 

0.8h 


0 

0.132 

0.238 

0.373 

0.526 

-0.012 

0.017 

0.065 

0.140 

0.5 ‘ 

0.135 

0.241 

0.376 

0.529 

-0.010 

0.019 

0.068 

0.143 

0.75 

0.137 

0.243 

0.377 

0.529 

-0.008 

0.021 

0.069 

0.145 

1.0 


0.243 

0.377 

0.529 


0.023 

0.071 

0.146 

1.5 


0.242 

0.374 

0.526 


0.027 

0.074 

0.148 

2.0 



0.367 

0.519 



0.075 

0.148 

4.0 



0.313 

0.459 



0.072 

0.132 

6.0 




0.386 




0.108 

8.0 




0.326 




0.088 




r 
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Table 4. Normalized stress intensity factor Ic/ICq at the deepest 
penetration point y=0, L=Lq of an inner semi -elliptic 
circumferential surface crack in a cylindrical shell under 
uniform membrane loading Nco. 





u.icn 


J 


‘‘0 

u.t+n 


^2 

mm 

a=2h 

a=4h 

a=8h 


msm 

a=2h 

a=4h 

a=8h 

0 



n 

0.961 


0.507 

0.627 

0.741 

0.837 

0.5 


0.879 


0.960 


0.497 

0.618 

0.735 

0.833 

0.75 

0.804 

0.875 


0.959 


0.487 

0.610 

0.729 

0.829 

1.0 


0.870 

0.923 

0.957 



0.600 

0.722 

0.824 

1.5 


0.858 

0.916 

0.953 



0.579 

0.704 

0.812 

2.0 



0.907 

0.948 




0.685 

0.798 

4.0 



0.870 

0.926 




0.613 

0.739 

6.0 




0.902 





0.687 

8.0 




0.881 





0.646 



'-0 = 

0.6h 




1-0 - 

0.8h 


0 

0.245 


0.451 

0.582 




0.149 

0.216 

0.5 

0.240 

0.330 

0.444 

0.576 




0.147 

0.213 

0.75 

0.236 

0.324 

0.438 

0.570 


0.073 


0.145 

0.210 

1.0 


0.318 

0.431 

0.563 



0.101 

0.143 

0.207 

1.5 


0.305 

0.414 

0.546 



0.101 

0.140 

0.200 

2.0 



0.398 

0.529 




0.137 

0.194 

4.0 



0.350 

0.467 




0.133 

0.177 

6.0 




0.422 





0.168 

8.0 




0.392 





0.163 
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Table 5. Nonnalized stress intensity factor k/kg at the deepest 
penetration point y=0, L=Lq of an inner semi -elliptic 
circumferential surface crack in a cylindrical shell 
under uniform bending moment Moo. 

I 




"o = 

u.zn 



*-0 = 

0.4h 


^2 

a=h 

a=2h 

a=4h 

a=8h 

a=h 

a=2h 

a=4h 

a=8h 

0 

0.804 

0.875 

0.926 

0.959 

0.441 

0.579 

0.710 

0.819 

0.5 

0.797 

0.870 

0.923 

0.957 

0.429 

0.569 

0.703 

0.814 

0.75 

0.789 

0.866 

0.921 

0.956 

0.418 

0.559 

0.696 

0.809 

1.0 


0.860 

0.917 

0.954 


0.547 

0.687 

0.803 

1.5 


0.847 

0.909 

0.950 


0.522 

0.666 

0.789 

2.0 



0.900 

0.945 



0.643 

0.772 

4.0 



0.859 

0.920 



0.557 

0.702 

6.0 




0.894 




0.640 

8.0 




0.871 




0.592 



‘-0 = 

0.6h 



l-o' 

0.8h 


0 

0.132 

0.238 

0.373 

0.526 

-0.012 

0.017 

0.065 

0.140 

0.5 

0.125 

0.230 

0.364 

0.518 

-0.013 

0.015 

0.062 

0.136 

0.75 

0.119 

0.222 

0.356 

0.511 

-0.013 

0.014 

0.060 

0.133 

1.0 


0.214 

0.347 

0.502 


0.013 

0.057 

0.129 

1.5 


0.197 

0.326 

0.481 


0.012 

0.053 

0.120 

2.0 



0.306 

0.460 



0.049 

0.112 

4.0 



0.244 

0.382 



0.042 

0.089 

6.0 




0.327 




0.078 

8.0 




0.289 




0.070 


r 
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Table 6. Normalized stress intensity factor k/ko at the deepest 
penetration point y=0, L=Lq of an outer semi-elliptic 
axial surface crack in a cylindrical shell under uniform 
membrane loading No, 




■-0 - 

u.^n 



"o - 

u.^fn 


^1 

msm 

a=2h 

a=4h 



a=8h 

msm 

a=2h 

a=4h 

a=8h 

0 

0.817 


B1 

Bi 

0.507 

mi 

0.741 

0.837 

0.5 

0.822 

|H 


mm 

0.518 


0.748 

0.841 

0.75 

0.826 

0.888 

IH 

iig 

0.527 

|H 

0.752 

0.844 

1.0 


0.890 

0.934 

0.963 


0.649 

0.757 

0.847 

1.5 


0.894 

0.936 

0.964 


0.663 

0.766 

0.853 

2.0 



0.938 

0.965 



0.773 

0.857 

4.0 



0.935 

0.964 



0.775 

0.860 

6.0 




0.959 




0.848 

8.0 




0.954 




0.834 



'-0-- 

0.6h 



Lo = 

0.8h 


0 



0.451 

ra 

0.073 


0.149 

0.216 

0.5 

0.255 

0.346 

0.461 

Hi 

0.076 


0.154 

0.223 

0.75 

0.264 

0.355 

0.468 

0.597 

0.080 


0.159 

0.229 

1.0 


0.364 

0.477 

0.604 


0.118 

0.165 

0.235 

1.5 


0.384 

0.494 

0.619 


0.130 

0.178 

0.250 

2.0 



0.509 

0.631 



0.192 

0.264 

4.0 



0.532 

0.651 



0.225 

0.299 

6.0 




0.641 




0.303 

8.0 




0.622 




0.294 
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Table 8. Normalized stress intensity factor k/ko at the deepest 
penetration point L=Lq> y=0 of an inner semi-elliptic 
axial surface crack in a cylindrical shell under uniform 
membrane loading N<» 





■-0 - 

u.cn 



■-0 - 



^1 

■IHI 


a=4h 

a=8h 

a=h 

a=2h 



a=4h 

a=8h 

0 

0.817 



0.961 

0.507 

B 

0.741 

0.837 

0.5 

0.813 

0.880 

mm 

0.960. 

0.501 

B 

0.737 

0.834 

0.75 

0.810 

0.878 

mm 

0.960 

0.498 

mm 

0.734 

0.832 

1.0 


0.876 

0.926 

0.959 


IM 

0.732 

0.830 

1.5 


0.873 

0.924 

0.958 


0.611 

0.728 

0.827 

2.0 



0.922 

0.957 



0.725 

0.825 

4.0 



0.916 

0.953 



0.718 

0.819 

6.0 




0.950 




0.811 

8.0 




0.946 




0.802 



•■0 = 

0.6h 



>•0 = 

0.8h 


0 



0.451 



0.104 

0.149 

0.216 

0.5 

0.243 

0.333 

0.447 

m 

0.074 

0.104 

0.148 

0.215 

0.75 

0.243 

0.331 

0.445 

0.576 

0.075 

0.105 

0.149 

0.215 

1.0 


0.331 

0.443 

0.574 


0.107 

0.150 

0.215 

1.5 


0.333 

0.444 

0.572 


0.112 

0.153 

0.217 

2.0 



0.444 

0.571 



0.158 

0.221 

4.0 



0.451 

0.570 



0.177 

0.237 

6.0 




0.569 




0.242 

8.0 




0.561 




0.241 
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Table 9. Normalized stress intensity factor k/ko at the deepest penetration 
point L = Lo> y = 0 of an inner semi -elliptic axial surface crack 
in a cylindrical shell under uniform bending moment Mo,. 




‘-o 

u.zn 



■-0 - 

U. 


^1 

a=h 

a=2h 

a=4h 

a=8h 

a=h 

a=2h 

a=4h 

a=8h 

0 

0.804 

0.875 

B 

0.959 


0.579 

0.710 

0.819 

0.5 

0.799 

0.872 

mm 

0.958 


0.573 

0.706 

0.815 

0.75 

0.796 

0.869 

0.923 

0.957 

0.430 

0.568 

0.702 

0.813 

1.0 


0.867 

0.921 

0.956 


0.565 

0.699 

0.811 

1.5 


0.864 

0.919 

0.955 


0.560 

0.694 

0.807 

2.0 



0.917 

0.954 



0.691 

0.805 

4.0 



0.911 

0.950 


. 

0.682 

0.797 

6.0 




0.946 




0.788 

8.0 




0.942 




0.777 



'-0 = 

0.6h 




0.8h 


0 

mi 

0.238 


0.526 


0.017 


0.140 

0.5 

|H 

0.233 

0.368 

0.521 


0.017 


0.139 

0.75 

0.128 

0.231 

0.365 

0.518 

-0.010 

0.018 

0.064 

0.138 

1.0 


0.230 

0.363 

0.516 


0.019 

0.065 

0.138 

1.5 


0.233 

0.363 

0.513 


0.024 

0.069 


2.0 



0.363 

0.513 



0.074 

0.145 

4.0 



0.369 

0.515 



0.091 

0.161 

6.0 




0.507 




0.164 

8.0 ■ 




0.495 




0.161 


- 24 - 




















Table 10. Distribution of the normalized stress intensity 

factor k/ko along the crack front in a cylindrical 
shell containing an inner or outer semi-elliptic 
circumferential surface crack (see insert in 
Fig. 4), X 2 = 2, a=4h, Lo=0.4h, v=0.3. 


7T 

Outer Crack i 

1 

Inner Crack 

Membrane 
Loadi nq 

Bending j 

Membrane 

Loading 

Bending 

1.0 

0.727 

0.692 

0.685 

0.643 

0.894 

0.719 

0.689 

0.678 

0.641 

0.789 

0.694 

0.680 

0.658 

0.637 

0.684 

0.655 

0.665 

0.625 

0.628 

0.578 

0.604 

0.643 

0.580 

0.615 

0.473 

0.544 

0.618 

0.527 

0.597 

0.367 

0.477 

0.583 

0.465 

0.569 

0.263 

0.406 

0.538 

0.399 

0.529 


Table 11. Distribution of the normalized stress intensity factor 
k/ko along the crack front in a cylindrical shell con- 
taining an inner or outer axial semi-elliptic surface 
crack (see insert in Fig. 4), v=0.3. 



Inner Crack 
a=h,R^- = 10h,LQ=0.2h 

Inner Crack 
a=4h,Ri=10h,Lp=0.8h 

Outer Crack 
a=4h,LQ=0.4h,X-|=2 

TT 

Tension 

Bending 

Tension 

Bending 

Tension 

Bending 

1.0 

m 

0.799 


0.078 

0.773 

0.747 

0.894 

msm 

0.797 



0.764 

0.743 

0.789 


0.792 

0.157 

0.094 

0.736 

0.731 

0.684 

mm 

0.782 


0.109 

0.693 

0.710 

0.578 

0.730 

0.765 

0.147 

0.124 

0.637 

0.683 

0.473 

0.685 

0.739 

0.138 

0.139 

0.572 

0.652 

0.367 

0.628 

0.700 

0.126 

0.149 

0.500 

0.612 

0.263 

0.559 

0.642 

0.114 

0.154 

0.426 

0.562 
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Figure 2. Comparison of the stress intensity factors obtained from 
the line-spring shell model and the axi symmetric elas- 
ticity solution [20]. (a) Stress intensity factor at 

the deepest penetration point of an external semi-elliptic 
circumferential crack in the shell, (b) same as (a) for 
an internal surface crack, (c) elasticity solution for the 
external axisymmetric crack, (d) the internal axisymmetric 
crack. (For Lg=0.6h, kQ=4.035 ao>^» kp=0.582 ko> 
oq: uniform axial stress, a=8h) 



4 
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Figure 3. Comparison of the line-spring shell stress intensity 
factor at the deepest penetration point of an internal 
axial surface crack (dashed lines) with the corresponding 
plane strain ring solution (full lines) [18j^a=8h 





Figure 4. 


Variation of the stress intensity factor along the front 
of a semi-el liptic external circumferential surface 


crack in a cylindrical shell. Xj, = 2, a .= 4h, Lo = 0.4h 
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Figure 5. 


Same as Figure 4, for internal surface crack (Ao= 2, 
a=4h, Lo=0.4h). 


Figure 6. Variation of the stress intensity factor for a semi- 
elliptic internal axial surface crack in a cylindrical 





Figure 7. Comparison of the line spring shell results (dashed lines) 
with the finite element solution (full lines) [14] for a 
pressurized cylinder containing a semi -elliptic internal 
axial surface crack. 
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